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Abstract 

Neutrices are additive groups of negligible functions that do not contain any constants except 
0. Their calculus was developed by van der Corput and Hadamard in connection with asymptotic 
series and divergent integrals. We apply neutrix calculus to quantum field theory, obtaining finite 
renormalizations in the loop calculations. For renormalizable quantum field theories, we recover 
all the usual physically observable results. One possible advantage of the neutrix framework is 
that effective field theories can be accommodated. Quantum gravity theories appear to be more 
manageable. 
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I. INTRODUCTION 



Quantum field theory is the proud product of quantum mechanics and special relativity. 
It is difficult to contemplate modern theoretical physics without it, for the very simple reason 
that it is extraordinarily successful. Coleman, a master of quantum field theory, compares 
the triumph of quantum field theory to "a glorious triumph parade, full of wonderful things" 
that make "the spectator gasp with awe and laugh with joy."[l| But there is a fly in the 
ointment: in general, infinities pop up in loop corrections in quantum field theory. In 
renormalizable theories, these infinities can be renormalized away in a process that Feynman 
likened to sweeping dirt under the rug. Even quantum electrodynamics, arguably the most 
phenomenologically successful theory in physics, has this "defect." Dirac was quoted to 
have said, "Quantum electrodynamics is almost right." He was bemoaning the divergent 
integrals that plague loop calculations in QED. And it is probably not wrong to think 
that Schwinger put forth his source theory approach to replace the conventional quantum 
field theory partly because of the ultraviolet divergences in the latter, To put it simply, 
quantum field theorists are confronted with the challenge to rid quantum field theory of 
infinities while preserving its many spectacular successes in the process. 

In another development, Dyson showed that the power series in the electron charge e for 
QED cannot be a convergent series. |2] His argument went roughly as follows. If the series is 
convergent, then such a series has a finite radius of convergence in the complex plane around 
e = 0. This means that one can replace e by ie without encountering a discontinuity. But the 
physics under this replacement is actually very discontinous: electrons would attract each 
other, likewise for positrons; on the other hand, opposite charges would repel each other so 
that electrons would move away as far as possible from positrons. Then every physical state 
would be unstable against the spontaneous creation of large numbers of particles. Thus the 
series for QED cannot be a convergent series. Dyson went on to suggest that the series is an 
asymptotic series. Taking his argument at face value, one should look for a proper tool to 
handle asymptotic series, for perturbative QED in particular and perturbation calculations 
in quantum field theory in general. 

Better yet, we should look for a mathematical tool that can handle asymptotic series in 
quantum field theory and rid it of infinities at the same time. WeQ] have suggested that 
such a tool is already available in the neutrix calculus developed by van der Corput[5j and 
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Hadamard. But instead of beginning with the axioms of neutrix calculus, in this paper we 
follow the path suggested to us by Prof. E. M. de Jager. In the next section, we discuss 
the close relation between asymptotic series and the finite parts of divergent integrals, and 
then introduce neutrix calculus to regularize divergent integrals. To make the presentation 
easy to follow, we use various examples to illustrate the approach. In Section III, we apply 
neutrix calculus to one-loop calculations in QED and 4 theories. In Section IV, switching 
to dimensional regularization, we apply neutrix calculus to the pure Yang-Mills system and 
"rederive" asymptotic freedom in one-loop. In both these two sections, there is actually 
very little new calculation involved, but they do serve to illustrate how neutrix calculus 
converts infinite renormalizations to mathematically well-defined finite renormalizations. 
As far as we can tell, all finite physically meaningful results of renormalizable quantum field 
theories are recovered in the neutrix framework. Section V is devoted to conclusions and 
a short discussion on the application of neutrix calculus to quantum gravity and effective 
quantum field theories in general. For completeness, we include two appendices to further 
discuss some mathematical properties of neutrices. In Appendix A, we consider more general 
Hadamard neutrices and a more general type of asymptotic series. To do quantum field 
theory in configuration space, one multiplies operator- valued distributions of quantum fields. 
In Appendix B, we show that the use of neutrix calculus allows one to put these products 
on a mathematically sound basis. 

II. FROM ASYMPTOTIC SERIES AND DIVERGENT INTEGRALS TO NEU- 
TRIX CALCULUS 

Neutrix calculus gives us a framework to deal with asymptotic series and divergent inte- 
grals. To motivate its introduction we follow de Jager's suggestion to begin with a standard 
example of asymptotic series. Consider 



for x — > oo. The "normal" procedure is to expand the denominator and write formally 




(1) 




(2) 
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The right hand side is clearly a divergent series for any finite value of x. This is not surprising 
because the expansion of the denominator is allowed only for t < x, but the integration over 
t extends to oo. 

Nevertheless one can distill from Eq. (JTJ a different type of expansion which avoids the 
rigid rules of the above power series approach. Here we find the asymptotic series by re- 
peatedly integrating by parts: 



e-* 



/(*) = 

x + 1 

giving 



oo e ~t 



o (x+ty 



T dt, (3) 



1 1 , (n — ~\V r°° <°~* 

f( x ) = _ + ... + Zl + (_i)« n ! / _ dt. (4) 

Jy 1 x x 2 y ' x n y ' Jo (x + t) n+1 y 1 

This equation is exact and is valid for any finite n. If we stop at the nth term, the remnant 
is given by 

(-!) n " ! I" , 6 \ o-i dt = °( -^tt) ■ (5) 
v ; Jo {x + t) n + l \x n+1 J v ; 

We get a satisfactory approximation for a sufficiently large x (x — n, e.g., would do). This 
is an asymptotic series, not a convergent series. If we stop at the nth term, we will have 
subtracted from the series the infinite part which, formally, is given by 1) J jx~^ +1 \ 

Note that in the series expansion the coefficients grow beyond bound for n — > oo. This is 
overcome by going to a larger x, a freedom not allowed for the "normal" (convergent) power 
series expansion in x. (The phenomenon of the coefficients growing with n happens often 
for asymptotic series, but it is not a definition of such a series.) 

Another example of asymptotic series is provided by the Stieltjes series J§] for the Stieltjes 
integral 

y(x) = / dt 

yy 1 Jo 1 + xt 

= l-x + 2!x 2 + ... + (-l) n n!x" + e n (x), (6) 



where 



and 



t ( x ) = {-l) n+1 {n + l)\x n+1 / (1 + xt)- n - 2 e- l dt, (7) 

Jo 



\e n (x)\<(n + l)\x n+1 «x n , x^0 + . (8) 

Next we move on to divergent series and their finite parts (i.e., Hadamard's "Parti Fini"). 
Here we encounter a similar situation like that treated above for asymptotic series. For this 



discussion, we use examples that will be relevant to the one-loop field theory calculations in 
the next section. Let us consider 

roc rp 

F' = / — — dx. (9) 

Jo 1 + x w 

Partitioning the integration into two parts, we can write 

pa rp rOO rp 

F' = / — — dx + / — — dx, (10) 

JO 1 + X J a 1 + X 

for any finite a. The second term contains the infinite part and the approximation of F' by 
the finite part is improved by increasing a. But we cannot take a = oo. We notice that here 
a corresponds to the n in Eq. (J3J). Doing the integration for the finite part, we obtain 

/•OO rn 

F' = a- log(l + a) + / dx, (11) 

J a 1 + X 

valid for all finite values of a. Increasing a indefinitely and subtracting the infinite part we 
get 

POO rp 

F' = / — — dx — > 0. (12) 
Jo l + x K J 

As another example, we consider 

x 



dx, (13) 



lo (l + x) 2 

which is also relevant to the physics we discuss in the next section. Proceeding as in the 
calculation of F' we have 

a dx r a dx f 00 x , 

+ / ~, z^dx 



o (l + x) Jo (l + x) 2 J a (l + x) 2 

1 /"CO 

loga + - 1+ / —dx. (14) 

1 + a A (l + x) 2 v ' 



Subtracting the infinite part with a — > oo yields 

x 



dx — > -1. (15) 



lo (l + x) 2 

For our purposes, the above two examples illustrate and encapsulate the essence of neutrix 
calculus developed by J. G. van der Corput (partly based on Hadamard's earlier work) in 
connection with asymptotic series and divergent integrals. By definition, a neutrix is a class 
of neglgible functions which satisfy the following two conditions: 
(1). A neutrix is an additive group; 



(2). It does not contain any constant except 0. 

lowing van de r Corpuf, lerffl, we illuStr ate tee eeocepts ef neutrix ea,cu luS with the 
help of various examples. 

Example 1. The neutrix consists of functions which go to as x — > oo. We call this 
neutrix e(x). 

Example 2. This concerns the class of polynomials in x with real coefficients. The neutrix 
is the subgoup of those polynomials divisible by x 2 + 1. The procedure gives the complex 
numbers (and this is apparently how Cauchy defined complex numbers). 

Eample 3. The neutrix consists of functions 

u(0 = U(0 + e(0, (16) 

where e(£) — > as £ — > oo as defined in Example 1, and 

f/(0 = ci£ + c 2 £ 2 + ... + c log£ + ... (17) 

This neutrix is denoted by H^. 

Example 4. The following example was considered by Hadamard: 

r2 f Cs- l 2 s - Cs- 1 ^ for s ^ 

C / x s ~ l dx = I (18) 
7 « I Clog2-Clog£ fors = 

Here £ > 0, s is a real number, and C is complex. For s > 0, the integral converges. For 
s < 0, s" 1 ^ grows without bound as £ — > 0, and so does — log£ for s — 0. Hadamard 
called Cs~ 1 2 s and CTog2 the finite part ("parti fini"), and Cs _1 £ s and Clog£ the infinite 
part ("parti infini"). He neglected the infinite part. For £ —>■ 0, van der Corput wrote the 
integral as 

/•2 1 f Cs-^ 5 for s 4 

C7 / x s - 1 dx= { (19) 
[ Clog2 for s = 0. 

Note the analytic extension in the complex s plane of the answer for Re s > to the entire 
complex plane with the exclusion of s = 0. 

In terms of the neutrix defined in Example 3, we can rewrite Eqs. ()12|) and (fT5Jl as 

/ -^—dx = 0, (20) 

Jo 1+x K ' 

and 

-dx = -1, (21) 



[1 + x) 



respectively. 

We conclude this section by pointing out that neutrices are essential to find the coefficients 
of an asymptotic series of the type 

f(x) = ... + c 2 x 2 + cix + c + c i \ogx + ... + c^ix' 1 + c_ 2 a;" 2 + ... + c_i/x _1 logx + ... (22) 

To find the coefficient c_i, for example, we multiply both sides by x and take the Hadamard 
Ho limit at x = 0. On the right hand side of Eq. (|22jl . the terms to the left of c_i are part 
of e(0) and hence converge to zero. On the other hand, the terms to the right of c_i are all 
in Hq and hence are neglected. We conclude 

c_i = limx/(x). (23) 

More mathematical properties of neutrices and asymptotic series can be found in the two 
appendices. 



III. NEUTRIX CALCULUS APPLIED TO 1-LOOP QED AND ^ THEORY 

As the first example in the application of neutrices to quantum field theory, let us consider 
QED: 

£ = ■ {-d - eA) + 77# - h,d,A u - d u A,) 2 . (24) 
We begin with the one-loop contribution to the electron's self energy in the Feynman gauge, 

=M = / 7f ^ h7 Jr',\: mh : . (25) 



(2tt) 4 [k 2 + \ 2 }[(p - k) 2 + m 2 ] 

where m is the electron bare mass and we have given the photon a fictitious mass A to 
regularize infrared divergences. Expanding S(p) about 7 -p — —m, 

E(p) = A + B(yp + m)+K, (26) 

one finds (cf. results found in Ref.Q) 

A=-7^rn(~D+~), (27) 



2tt V2 4 



a 
Air 



( f l dx 








\ m 





(28) 
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where a = e 2 /4ir is the fine structure constant and D is given by 

J_ r d 4 k 



in 2 J (k 2 + m 2 ) 2 
00 k 2 c\k 2 

(29) 



/o (A; 2 + m 2 ) 2 

with the second expression of D obtained after a Wick rotation. We note that 71, the last 
piece of S(p) in Eq. (|26|). is finite. Mass renormalization and wavefunction renormalization 
are given by m ren = m — A and ^/v en = Z 2 respectively with Z^ 1 = 1 — 13. Now, 
introducing dimensionless variable q = k 2 /m 2 , we bring in the neutrix TTqo to write D as 

"-r^" 1 - <3o) 

where, for the last step, we have used Eq. (|2ip. Since D = — 1 is finite, it is abundantly 
clear that the renormalizations are finite in the framework of neutrix calculus. There is no 
need for a separate discussion of the electron vertex function renormalization constant Z\ 
due to the Ward identity Z\ = Z 2 - 

The one-loop contribution to vacuum polarization is given by 

n '- w = !e2 / wf r (vfe+V^Vfr-VJ ' (31) 

A standard calculation [t| shows that IL^ takes on the form 

n M „ = 5m \ v + (k\„ - k^)U(k 2 ), (32) 
where ^ is the flat metric (+++-), 

5m 2 = ^-(m 2 D + D'), (33) 

27T 

and 

rt ^ 9ft fl Ik 2 \ 

K(k 2 ) = ~{D + ~) + — / - a:) log 1 + ^rx(l - x) , (34) 

67T 7T JO \ 7TT / 

with 

and 73 given by Eq. ()29j1 . Just as 73 is rendered finite upon invoking neutrix calculus (see 
Eq. PHI), so is D'\ 

2 jrfg 
Q 



D> = m 2 "A = 0, (36) 



S 



where, for the last step, we have used Eq. (J20J). Thus neutrix calculus yields a finite renor- 

— 1/2 

malization for both the photon mass and the photon wavefunction A^ en = Z 3 (and 
consequently also for charge e ren = Z\j 2 e) where Z^ 1 = 1 — 11(0). In electron-electron scat- 
tering by the exchange of a photon with energy-momentum k, vacuum polarization effects 
effectively replace e 2 by e 2 /(l — II(/c 2 )), i.e., 

2 



2 , „2 



e — ► e 



' " 1 11(7.-' 

„2 



z 3 (i-n(P)) 

„2 



i - (u(k 2 ) -n(o))' 

Eq. (pMjl can be used, for k 2 ^> m 2 , to show that 



(37) 



= \ ^—f P V- ( 38 ) 

1 - ^ log — a 



Thus we have obtained the correct running of the coupling S] with energy- momentum in 
the framework of neutrices. In fact, the only effect of neutrix calculus, when applied to 
QED (and other renormalizable theories), is to convert infinite renormalizations (obtained 
without using neutrix calculus) to mathematically well-defined finite renormalizations. As 
far as we can tell, all (finite) physically observable results of QED are recovered. In passing 
we mention that the use of neutralized integrals does not affect the results of axial triangle 
anomalies. 

Let us now apply neutrix calculus to the </> 4 theory, 

£=-I^^-ImV- J^- (39) 

The one-loop self-energy is given by 

, A f d 4 k 1 . . 

^ = '2 J + (4 °» 

which, by Eqs. (135(1 and (|36(h vanishes! Thus there is no mass renormalization (m 2 en = 

m 2 — E(0) = m 2 ) and no wave-function renormalization (wave-function renormalization 

constant = 1) in the one-loop approximation. The four-point function T(s,t,u), for 

incoming particle momenta pi, receives the following one-loop contributions 

A 2 f d 4 p 1 1 

T(s,t,u) = — / -, — —r— r- + (2 crossed terms). (41) 

V '' ; 2 J (2tt)V + ™ 2 {p-q) 2 + m 2 K ) \ I 



Here q = pi + P2, and we have used the Mandelstam variables s = —q 2 , t = — (pi + -P3) 2 , 
and u = — (pi +P4) 2 . A straight-forward calculation yields 

r(>, t, u) = 1— {3D - J Q dz (log(l - ^[1 - *]) + (s - t) + (s - u)) } . (42) 

Thus, the coupling receives a finite renormalization 

3DA\ 



Aren — A ( 1 



Mi + ^-l , (43) 



327T 2 

3A 

where we have defined X ren by T(s, t, u) at s = t — u — 0, and have used D = —1, given by 
Eq. (|3*0]h for the last step. 



IV. DIMENSIONAL REGULARIZATION AND PURE YANG-MILLS THEORY 

As shown by the appearance of photon mass in the above discussion of vacuum polar- 
ization in QED, the application of neutrix calculus to the energy-momentum cutoff reg- 
ularization scheme is not too convenient to use for more complicated theories like those 
involving Yang-Mills fields. For those theories, one should use other regularization schemes 
that manifestly preserve the Ward identity. In this regard, we note that already in 1961, van 
der Corput suggested that, instead of finding the appropriate neutrices, one can continue 
analytically in any variable (presumably including the dimension of integrations) contained 
in the problem of tackling apparent divergences to calculate the coefficients of the corre- 
sponding asymptotic series. It so happened that this was the approach taken by 't Hooft 
and Veltman who spearheaded the use of dimensional regularizationsiS)]. Let us now explore 
using neutrix calculus in conjunction with the dimensional regularization scheme. In that 
case, negligible functions will include 1/e where e = 4 — n is the deviation of spacetime 
dimensions from 4. In the calculations, the internal energy-momentum integration is now 
over n dimensions. The one-loop contributions to the electron's self-energy in QED is given 
by 

nv) = ^r(2 - -) t dx (n-2) 1 -p(l-x)+nm 

yP> (4tt)"/ 2 1 2' Jo [p 2 x(l-x) + m 2 x + X 2 (l-x)] 2 - n / 2 ' 1 } 

Again we expand S(p) about 7 ■ p = —m as in Eq. (}2l)|) . Using the approximation for the 

gamma function 

Bmr( C ) = i-7 + 0(e), (45) 
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where 7 ~ .577 is the Euler-Mascheroni constant, and the approximation 

r~l + elog/, (46) 

for e <C 1, one finds 

A = —[3(7 - log47r) + 1] + — m / dx(l + x) logD , 
An In jo 

where Do = ^ 2 ^ 2 + A 2 (l — x), 1 and we have invoked neutrix calculus in dropping 
—3am/(2ne) and — a/(2ne) from A and i3 respectively. (In passing we note that TZ in 
the expansion of S(p) in the dimensional regularization scheme is the same as in the energy- 
momentum regularization used in the last section.) 

The one-loop vacuum polarization takes on the same form as given by Eq. (|32p. but now 
with 5m 2 = 0, and 

2 r ( 2 -D f 1 J_ s(l-x) 



U(k 2 ) = -8e 2 , K r . 2 / / rfx-r-- ^ hr- jr. (48) 

v ; (4tt)™/ 2 Jo [ft(l-s) + m 2 ]W2 v ; 



Ol 2ot 

Il(fc 2 ) = — [7 — log47r] H / dxxil — x) log[m 2 + x(l — x)k 2 

on n Jo 



In the limit e = 4 — n^Owe obtain 

2a r 1 
n Jo 

^ = l-^[7-log47r]--^logm 2 , (49) 
Z3 on on 

where we have followed neutrix calculus in dropping — e 2 /(6n 2 e) from II(A; 2 ). By design, 
the generalized neutrix calculus renders all the renormalizations finite. Again, all physically 
measurable results of QED appear to be recovered. 

Next we consider pure non-Abelian gauge theroy for a group G with structure constants 

C a bc 

C = -\{d»A v a - d v A» + g C abc A%A») 2 . (50) 



1 Since Dq is dimensionful, the appearance of logZ?o is very strange. This is due to the fact that the e 2 
in Eq. I|44(l is in fact also dimensionful for e 7^ 0. If we replace e 2 there by e 2 /i 4_n with a certain mass 
/1 so that the e 2 will then be dimensionless, Dq will be replaced by dimensionless Dq/h 2 . For simplicity 
we have not written this dimensional dependence out explicitly Hereafter this dependence will be left 
implicit. 



11 



The one-loop contribution to the self-energy of the gauge field, in the Feynman gauge, is 
given by 

n£(*) = (*V- Wn^k), (51) 

with 

U cd (k) = C 2 cd ^^T(2 - f Q dxV^[(A - 2n) + (6n - 4)s + (8 - 4n)x 2 ], (52) 

where the symbol C 2 cd stands for C a b c Cabd and V = k 2 x(l — x). 
With the aid of Eqs. (gnj) and (gHJ), we get 

n- W = ^ (fi - " log P + H,_ 7 + l 0g 47r] + 7-i). (53) 

Invoking neutrix calculus to remove the negligible e _1 term and specializing to the SU(N) 
group where C 2 d = NS c d, we obtain 

U cd (k)=5 cd U(k 2 ), (54) 

with 

n( ^M(_^ + ...), m 

where ... stands for the finite /c-indepentdent terms in Eq. (J53j) . For QED we have used 
the on-shell renormalizations with massive electrons. For pure Yang-Mills theory we will 
do the renormalizations for massless gauge fields at space-like energy-momentum /i 2 . The 
renormalization constant Z 3 is then given by 

Z3= i^TTp)" 1_ TToV log/i + - (56) 

Similarly the renormalization constant for the three-point function in the Feynman gauge 
can be calculated to yield 



Z 1 ^l--^log/* + ... (57) 



4iV^ 

3167T 

where ... stands for finite /i-independent terms. The Callan-Symanzik function for the renor 

1 3/9 

malized coupling g — Z 1 Z 3 g in the one-loop approximation is then given by 



11 No 3 . . 
- r (58) 



3 16tt 

Thus we have recovered the negative f3 function for SU(N) Yang-Mills theory. 
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V. CONCLUSIONS AND DISCUSSIONS 



In this paper, we have proposed to apply neutrix calculus, in conjunction with Hadamard 

rn 

integrals, developed by J.G. van der Corput [5| in connection with asymptotic series, to 
quantum field theory, to obtain finite results for the coefficients in perturbation series. The 
replacement of regular integrals by Hadamard integrals in quantum field theory appears to 
make good mathematical sense, as van der Corput observed that Hadamard integrals are 
the proper tool to calculate the coefficients of an asymptotic series. (Actually Hadamard 
integrals work equally well for convergent series.) 

For renormalizable quantum field theories like QED, 4 theory, and pure Yang-Mills, we 
have demonstrated (even if not adequately) that we recover all the usual physically observ- 
able results. The only effect of neutrix calculus appears to just change the "amount" of 
renormalizations — from infinite renormalizatins to finite renormalizations. But what about 
non-renormalizable field theories like quantum gravity? Using dimensional regularization, 
'tHooft and Veltman[9] found that pure gravity is one-loop renormalizable, but in the pres- 
ence of a scalar field, renormalization was lost. For the latter case, they found that the 
counterterm evaluated on the mass shell is given by ~ e~ l ^/~gR 2 with R being the Ricci 
scalar. Similar results for the cases of Maxwell fields and Dirac fields etc (supplementing 
the Einstien field) were obtained^]. It is natural to inquire whether the application of 
neutrix calculus could improve the situation. Our preliminary result is that now essentially 
the divergent e _1 factor is replaced by —7+ finite terms. 

Since neutrix calculus does not tolerate infinities, we conjecture that it can also be used 
to ameliorate other problems that can be traced to ultraviolet divergences. Two well-known 
problems come to mind. The hierarchy problem in particle physics is due to the fact that 
the Higgs scalar self-energies diverge quadratically. It is conceivable that neutrix calculus 
may help. The cosmological constant problem can be traced to the quartic divergences in 
zero-point fluctuations from all quantum fields. Perhaps neutrix calculus can help with this 
problem too. Indeed, for a theory of gravitation with a cosmological constant term, the 
cosmological constant receives at most a finite renormalization from the quantum loops in 
the framework of neutrix calculus. 

By ridding quantum field theories of their ultraviolet divergences, neutrix calculus may 
appear to have fallen victim to its own successes, for now we have lost renormalizability as a 
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physical restrictive criterion in the choice of sensible theories. However, we believe that this 
is actually not as big a loss as it may first appear. Quite likely, all realistic theories now in our 
possession are actually effective field theories. 0, They appear to be renormalizable field 
theories because, at energies now accessible, or more correctly, at sufficiently low energies, all 
the non-renormalizable interactions are highly suppressed. By tolerating non-renormalizable 
theoriess, neutrix calculus has provided us with a more flexible framework to study all kinds 
of particle interactions. 

To the extent that it is relevant for asymptotic series and lessens the divergences in 
quantum field theories, neutrix calculus is a powerful tool that is, in our opinion, too valuable 
not to make use of. How much it can really help to ameliorate problems related to ultraviolet 
divergences in quantum field theory remains to be seen. But based on our study so far, we 
tentatively argue that neutrix calculus has banished infinities from quantum field theory. 
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Appendix A: The Hadamard neutrix H a and Hadamard series expansion 

In Section II, we considered some examples of asymtotic series. Essentially, the series 
f(x) = a + a,i(x — b) + 02 (x — b) 2 + ... for finite b is an asymptotic series [f| if and only if 
there exists an no > 0, such that for n > uq, 

Km - — ^— |/(a;) - a Q - a t {x - b) - . . . - a n (x - b) n \ = 0, (59) 

x^b [X — 0) n 
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with the remnant being at most ~ {x — b) n+1 . 

In this appendix, we consider the Hadamard series expansion and a generalized definition 
of asymptotic series. Let k be a region in the complex plane with a finite limit point a which 
does not belong to k. A crucial property is that a point £ in k can approach a in such a 
way that arg(£ — a) has a finite limit. We say that U(£) in k has a Hadamard expansion in 
powers of £ — a if £/(£) is defined in k, and for points £ near a, has an asymptotic expansion 
of the kind 

U(0 ~ E£** fc (£ - a)*- log fc H£ - a). (60) 

Here x/i, ^/i and integers k h >0 are independent of £, Re ^ — > oo as /i — >■ oo, and log fch (£) 
stands for (\og(£)) kh . 

The definition of asymptotic series here is more powerful than the one we used in Section II 
in that, for every number q (independent of £), there exists an integer mo > with the 
property that, for every number m > m , there are two positive numbers c m and e m such 
that, for every point £ in k with |£ — a\ < e m , 

U(0 - K=oXh (£ - af- log fc ^ - a) | < c m |£ - a|* (61) 

The new definition of asymptotic series is essential in deriving the analytic properties of 
functions represented by neutralized integrals. Note that this means that there is a radius 
of convergence which goes to zero as m goes to infinity. If, in Eq. ([61 J), there is no term with 

= kh = 0, then we speak of a Hadamard series expansion without constant term. 

We now define the Hadamard neutrix H a . It consists of all functions which in the 
neighborhood of a can be written as 

u(t) = U(£) + e(t). (62) 

Here U(£) has a Hadamard series expansion in powers of (£ — a) without a constant term 
and e(£) approaches zero when ^ in k approaches a. This neutrix is called H a ; it has a 
domain a variable £, and a is called the carrier. Thus every negligible functions in H a 
can be written as 

V&oXh (£ - log fch (e - a) + e(fl, (63) 

with no = kh = term. 

We end Appendix A with three examples. 
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Example 1. If two different points a and b are connected by a rectifiable curve in the 
complex plane, which has a tangent in a but does not contain a, then for every integer k > 
and for every complex s 

rb . , f i l g fc+1 (6-a) for s = 

Example 2. For every complex s, except 0, —1, —2, 

x s - 1 e~ x dx = r(s). (65) 



Ho 

Example 3. For — 7r < arg I < n, s not an integer < 0, and s + 1 not an integer > 0, one 
has 



fe 

///, 1 ' A fc+l v^.sj n /) 

Here the integral is along the positive axis. This integral contains the special case (with 
A = 2, k = 0, t = -a, and 2s = (3 + 1) 

J ho {x 2 + iy dx ~2 r^)FPP ' (67) 

an expression encountered in dimensional regularization calculations for which (3 = 3 — e 
with a positive infinitesimally small e. 



Appendix B: Neutrices and products of singular functions 

It has often been claimed that infinities find their way into quantum field theory because 
one works with products of singular functions (like Feynman propagators) in the theory. In 
what used to be called axiomatic field theory, one had to search for the definitions of the 
products of operator-valued distributions. But whereas generalized functions usually cannot 
be multiplied in the theory of distributions developed by Schwartz, the Hadamard approach 
adopted in this paper allows the multiplication for a wide class of distributions as we will 
show in this Appendix. In fact this is the reason that one obtains finite rather than infinite 
renormalizations in the framework of neutrix calculus. 

Lighthill^J has already defined and listed the products for some distributions in his 
book. Van der Corput 5] has been able to simplify and extend Lighthill's list by doing 
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what one usually does in quantum field theory, namely, by using Fourier transformed (FT) 
functions 

/oo 
e- tpx f(x) dx. (68) 
-co 

Then he defines the FT of the product of two functions by folding (as usual), but using 
Hadamard integtrals (instead of the "normal" integrals): 

FT[h{x) f 2 {x)\ = [ H °° ^F,{k) F 2 {p - k). (69) 

Consider, for example, the one- dimensional Dirac delta function multiplying itself 5(x) x 
S(x). In Schwartz' approach, this product is not mathematically meaningful because its 
Fourier transform diverges: 

/°° dk 
^lxl— oo, (70) 
-co ZTT 

where we have used the convolution rule in Fourier transform and have noted that the 
Fourier transform of the Dirac delta function is 1. In contrast, the Hadamard- neutralized 
Fourier transform of the product 5(x) x 5(x), 



L 



d,k 



yields 

5(x) x 5(x) = 0, (72) 

a mathematically meaningful (though somewhat counter-intuitive) result! One can also show 
that (for non-negative integers a and b) 

{d a 5{x)) (d b 5{x)) = 0, (73) 

because 

/ k a (p-k) b dk = 0. (74) 

Similarly one can show that, for non- negative integers m and a, the Fourier transform of 
x~ m d a 5(x) vanishes, and the Fourier transform of x m d a 5(x) is zero for m > a and i m+a a(a — 
I)... (a — m + \)p a - m for m < a respectively. 

Lastly let us generalize the above discussion for products of two 1-dimensional Dirac delta 
functions to the case of (3 + l)-dimensional Feynman propagators 



(27r) 4 p 2 + m 2 — ie 
= i-5(x 2 ) - —^=Hi 2 \mV=^ie), (75) 
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where is the Hankel function of the second kind and we use the (+++-) metric. The 
Fourier transform of A + (x) x A + (x) (which appears in certain quantum loop calculations, 
e.g., in Eq. (|42|) in the text) is given by 

f { d A k 1 1 

J d A xe ipx A + (x)A + (x) = J (2i)4p + m2 _ j£ ( p _ k y + m 2_ ie 

= / „ D- , 1 r= f 1 dz log ( l + lL z (l-z) \ , (76) 
4(2tt) 2 4(2tt) 2 Jo \ m 2 v ') ' v ; 

where -D is given by Eq. (|29p. But -D is logarithmically divergent. Hence A + (x) x A + (x) 
is not mathematically well defined. On the other hand, this product is mathematically 
meaningful in the Hadamard-van der Corput approach. To wit, D = — 1 in the neutralized 
version, hence A + (x) x A + (a?) ~ 5^(x)+ regular part, where S^'(x) is the 4-dimensional 
Dirac delta function. 



[1] S. Coleman, Aspects of Symmetry (Cambridge Univeristy Press, Cambridge, 1985). 

[2] J. Schwinger, Particles, Sources, and Fields, Volumes I, II, III (Addison- Wesley, Reading, 

Massachusetts, 1970, 1973, 1989). 
[3] F. J. Dyson, Phys. Rev. 85, 631 (1952). 
[4] Y. J. Ng and H. van Dam, |arXiv:hep-th/0410285 



[5] J. G. van der Corput, J. Analyse Math. 7, 281 (1959); MRC Technical Reports # 128 - 130, 
University of Wisconsin at Madison (1960); 1961 Amsterdam lectures (Voordrachten over 
neutrixrekening, Stichting Mathematisch Centrum, Amsterdam, 1961). For more references, 
see MathSciNet under "neutrix". 

[6] See, e.g., C. M. Bender and S. A. Orszag, Advanced Mathematical Methods for Scientists and 
Engineers (McGraw-Hill, New York, 1978). 

[7] J. M. Jauch and F. Rohrlich, The Theory of Photons and Electrons (Addison- Wesley, Cam- 
bridge, Massachusetts, 1955). 

[8] M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory (Addison- 
Wesley, Reading, Massachusetts, 1995); T. P. Cheng and L. F. Li, Gauge Theory of Elementary 
Particle Physics (Oxford University Press, Oxford, 1984). 



18 



[9] G. 't Hooft, Nucl Phys. B62, 44 (1973); M.J.G. Veltman in Quantum Theory of Gravitation, 
Les Houches session XXVIII (North Holland, Amsterdam, 1976). 

[10] See, e.g., S. Deser and P. van Nieuwenhuizen, Phys. Rev. Lett. 32, 245 (1974), and S. Deser, 
H.S. Tsao and P. van Nieuwenhuizen, Phys. Lett. 50B, 491 (1974). 

[11] S. Weinberg, The Quantum Theory of Fields, Volume I (Cambridge University Press, 1995). 

[12] M. J. Lighthill, Introduction to Fourier Analysis and Generalized Functions (Cambridge Uni- 
versity Press, Cambridge, 1959). 



19 



